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Abstract 

We obtain the solution of the Dirac equation in (2+1) dimensions in the presence of a con¬ 
stant magnetic field normal to the plane together with a two-dimensional Dirac-oscillator potential 
coupling. We study the energy spectrum of graphene quantum dot (QD) dehned by electrostatic 
gates. We give discussions of our results based on different physical settings, whether the cyclotron 
frequency is similar or larger/smaller compared to the oscillator frequency. This defines an effective 
magnetic field that produces the effective quantized Landau levels. We study analytically such field 
in gate-tunable graphene QD and show that our structure allow us to control the valley degener¬ 
acy. Finally, we compare our results with already published work and also discuss the possible 
applications of such QD. 
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1 Introduction 


In recent years, several interest have been devoted to the study of two-dimensional (2D) system such 
as quantum wells, quantum wires, and quantum dots [1-0]. This interest is due to the technological 
advances in nanofabrication. In addition, one of the most important recent development in semicon¬ 
ductor has been the achievement of structures in which the electronic behavior is essentially 2D. This 
means that the charge carriers are confined in a potential such that their motion in one direction is 
restricted and thus is quantized, leaving only a two-dimensional momentum. In particular, there has 
been considerable amount of work in recent years on semiconductor conhned structures, which finds 
applications in electronic and optoelectronic devices. The application of a magnetic field perpendicular 
to the heterostructure plane quantizes the energy levels in the plane, drastically affecting the density 
of states giving rise to the famous quantum Hall effect [7]. The latter remains as the most interesting 
phenomenon observed in physics because of its link to different theories and subjects. 

Graphene [8,9], two-dimensional crystalline materials, has became one of the most important 
subjects in condensed matter research in the last few years. This new material has a number of 
unique properties, which makes it one of the most promising materials for future nanoelectronics [10]. 
One of them is the band structure, which is gapless and exhibits a linear dispersion relation at two 
inequivalent points K and K’ in the vicinity of the Fermi energy. Moreover, its low energy of electrons 
are governed by a (2-1-1) dimensional Dirac equation. Those electrons behave as massless chiral 
fermions, i.e. relativistic electrons. Consequently, the electrons cannot be localized by any conhnement 
potential, which is related to the fact that electrons in graphene can have both positive and negative 
energies, i.e. Klein tunneling effect [11]. Graphene quantum dot (QD) [12-14] have sparked intense 
research activities related to quantum information storage and processing using the spin information 
of confined electrons. Various methods were used to make QD one of the most widely-used techniques 
using electrostatic gates [13]. 

On the other hand, the Dirac oscillator was proposed by Moshinsky and Szczepaniak [15] in 
1989 and is considered as the relativistic version of the harmonic oscillator. The Dirac oscillator 
has been studied extensively [16-21] because of their probable applications in many branches of 
physics. Additionally, the Dirac oscillator has been used in optics [22], Jaynes Gummings model [23] 
and grephene [24]. It is only recently that the first experimental microwave realization of the one¬ 
dimensional Dirac oscillator was developed [25]. 

We combine different approaches to achieve our goal. Indeed, based on [26, 27] we set the Hamil¬ 
tonian system of our problem where a similarity transformation is used to simplify the process for 
obtaining the solutions. Later on, we define the QD by gates introducing an electrostatic confining 
potential. We find the bound state solution of gate-tunable graphene QD in the presence of a constant 
magnetic field B and Dirac oscillator of frequency uj as well as a mass term that might be intro¬ 
duced by the underlying substrate [28,29]. The eigenspinors are obtained in terms of the confluent 
hypergeometric functions showing dependence of B and the oscillator coupling oj. 

Subsequently, we analyze the impact of the external field B on the solutions of the energy spectrum 
of the QD by extracting interesting properties. More precisely, we consider three different cases 
corresponding to the relative strength of B with respect to ui. In doing so, we start by defining 
an effective magnetic field, that produces the effective quantized Landau levels, and focus on its 
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dependence of the bound states in circularly symmetric QD. We show how to control the valley 
degeneracy by manipulating the effective magnetic field. This can help to form the valley filters, 
valves [30], or qubits [31], and spin qubits [32] in graphene. 

The paper is organized as follows. In section 2, we set our problem by reviewing some mathematical 
tools need to deal with our issues. To investigate the basic feature of the gate-tunable graphene QD, 
we set the appropriate confining potential and give the corresponding solutions of the energy spectrum 
in section 3. Using the matching condition at the boundary, we obtain the condition that governs the 
existence of the bound state. This will serve to study different liming cases related to the strength of 
the magnetic field. We conclude our results in the final section. 


2 Theoretical model 

In order to deal with our task we establish an appropriate Dirac equation describing our system. To 
go deeply in our study for the graphene QD, we introduce a mass term to open a gap. 


2.1 Dirac equation 

To start let us set some mathematical background related to Dirac formalism needed to deal with our 
task. Indeed, a particle of mass m in the presence of a constant perpendicular magnetic field can be 
described by consider the Dirac equation in (2 -|- 1) dimensions 

[z 7 ''(cl^ + - m ]-0 = 0 , // = 0 , 1,2 ( 1 ) 

where the electromagnetic potential = (^o,^) and the space-time gradient = (|,v). Here 
we have the representation 7 *^ = era and 7 = f cr with the 2x2 hermitian Pauli spin matrices {crj}^=i- 
The Dirac matrices 7 ^ satisfy the algebra 


7.} = = 0,1, 2 


( 2 ) 


with the metric Q = diag(-|- — —). 

In our study of system made of graphene, we need to consider massless Dirac fermions. To this 
end, we multiply (I) by (T 3 to open a gap, such as 



+ "mfi ] ip 


(3) 


where we have set a = ia^a and (3 = For time independent potentials, the two-component spinor 
wavefunction is separable 'p{t,r,9) = For regular solutions of (3), square integrability 

and the boundary conditions require that 'p{r,9) satisfies 


^/rp{r,9)\r=Q =0, ^{9 + 27^) = p{9). (4) 

r^oo 

To simplify the construction of the solution, we look for a local 2 x 2 similarity transformation A(r, 9) 
that maps the cylindrical projection of the Pauli matrices (cr • r, a ■ 9) into their canonical Cartesian 
representation (iti,cj 2 ), respectively. That is 

Aa • f = (Ti, KB ■ 9= CJ2. (5) 
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We note that any other choice for the pair of Pauli matrices can be obtained from the present one 
through a unitary transformation, hence leaving the physics of the problem unaltered. A 2 x 2 matrix 
that is defined by [26,27] 

A(r,6/) = 

We are interested to the Dirac oscillator for its probable application in many branches of physics as 
we noticed before [22,24,25]. Motivated by these investigations, we consider such oscillator in another 
context and emphasize its influence on a system based on the QD. To achieve this goal, we introduce 
an additional coupling as the 2D Dirac-oscillator potential [15,33], that keeps symmetry of the system. 
This coupling is introduced by the substitution V —)■ V + \ujrj3 where uj is the oscillator frequency 
and A is a constant parameter. To simplify the forthcoming analysis, we require that the condition 
X = m should be fulfilled. Now from the above consideration, we obtain the (2 + l)-dimensional Dirac 
equation for a charged spinor in static electromagnetic potential 

{H-E)x± = 0 (7) 


where the Hamiltonian is given by 


H = 


dr + iAr - Xujr - \dQ + Aq 


-dr - iAr - Xujr - ^dg + Ag 


+ Ao-3 + AqI 


( 8 ) 


and x± are the components of the transformed wavefunction \x) = A |V’), with A is given in (6). It is 
clearly seen that the second term is gap and third one can be regarded as an external potential. In 
the forthcoming analysis, we will fix different potential in order to deal with some basic features some 
properties of the gate-tunable graphene QD. 


2.2 Mass term 

To study the impact of the external held B and oscillating frequency a; we consider a system made of 
graphene described by the Hamiltonian (8) except we replace Ads by rAds with r = ±1 differentiates 
the two valleys K and K'. Thus, we call the wave function x± spinor. To go further, let us set 
some quantities such a constant magnetic held of strength B applied perpendicular to the (r, 0)-plane, 
which is B = B z. We choose the gauge A{r, 9) = ^ (0, r) and assume a circular symmetry in the 
conhnement potential Aq = U{r). 

Consequently, (7) becomes completely separable in radial and angular parts. Then, we can write 
the spinor wavefunction as 

X±ir,9) = (9) 

such that the angular component satishes the eigenvalue equation —= C T’ where ^ is a real 
separation constant giving the function 

r{e) = ( 10 ) 

V 27r 

On the other hand, the boundary condition 'ijj{9 + 2 tt) = 'ip{9) requires the following condition 

= +1 ( 11 ) 
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which, in turn, demand that = — 1 giving the quantum number 

e = k = ±1,±3,±5--- . (12) 

While, the Dirac equation for the two-component radial spinor is reduced to 

/ r\ + U-E ^ + £ + Gr\ / '/>+W 
V + ^ + Gr -tX + U - E j ^ 

where the physical constant G = X {ojc — oj) and the cyclotron frequency Thus, the presence of 

the 2D Dirac-oscillator coupling did, in fact, maintain the symmetry of the problem. It is interesting 
to note that the effect of the oscillator interaction is to produce an effective magnetic field. This 
effective magnetic field is assumed to originate in an effective form from the motion of the charge 
carriers relative to the planar hexagonal arrangement of carbon atoms. Moreover, its introduction is 
equivalent to change the magnetic field as 



B 


B - 2Xuj. 


(14) 


As a result of the wave equation (13), the two spinor components satisfy the relation 

1 


(fE {r) = 


E-U ±tX 


dir. 

=F:r + - + Gr 
dr r 




(15) 


Using (15) to eliminate one component in terms of the other in (13), which results in the following 
Schrodinger-like differential equation for each spinor component 

-{E-Uf + G (2^ ± 1)] I ctE^ir) = 0. (16) 

It is clearly seen that this equation can completely be solved by choosing an approriate form of 
the potential U{r). This will be done by specifying the nature of our system, which be will be the 
graphene QD one. However, for a constant potential and on the light of the above choice, we show 
that he solutions of (16) take the form 


and the function r(r) is given by the following combination 


r(r) = a^U 



a 

2 


+ l,\G\r^) 


+ PaM 



a 

2 


+ l,\G\r^) 


(18) 


with U and M are the confluent hypergeometric functions. The parameters Oo- and Pa are normal¬ 
ization constants, which can be fixed by the boundary conditions. We show that the corresponding 
eigenvalues e = E — U read as 


ria 


A2 - e2 1 

4\G\ ^2 




(19) 


where Ua is an integer value and s = sgnG. We notice that teh above results are in agreement with 
those reported in [26] for cj = 1. 
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3 Graphene quantum dot 


We show that how the above results can be used to study some physical properties of a system made 
of the graphene QD. This latter is defined by gates introducing an electrostatic confining potential 
as shown in Figure 1. In what follow we study the corresponding solutions in the presence of a 
constant mass term A, that will account for a gap 2A in the energy spectrum. In fact, we will consider 
three limiting cases those related to the strength of the external magnetic field with respect to Dirac 
oscillating frequency. 


3.1 Potential configuration 


We are interested by studying the impact of the Dirac oscillator coupling together with an external 
magnetic field on the graphene QD. In doing so, we look first for the bound states solution for the QD 
governed by the second order differential equation (16). To this end, we take a radially symmetric 
model with the following potential 


U{r) 


Uo, r > R 
0, r < R 


( 20 ) 


and its profile, which consists of two regions, is schematically shown in Figure 1. The radially sym¬ 
metric choice makes the analytic determination of the bound state solution possible. 


U(r) 


m 


▲ 







Uo 
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-► 


R 


Figure 1: Potential landscape of a QD with radius R. 


According to Figure 1, we distinguish between two regions in order to explicitly determine the solu¬ 
tions of the energy spectrum. Consequently, the general solutions to the radial equation (16), which are 
regular at the origin and decay exponentially as r —)■ oo, are given in terms of M (n^, |^ — f | -|- 1, |G| r^) 
inside the QD and U (n^, ~ f | + 1) 10*1 outside the QD. These are 


^ 2 (i+l«-^/ 2 |)/ 2 g-|G|rV 2 ^(|Wf 1+1/2) / “ f | + 1, |G| r^) , r> R 

\ /3.M(n<,|e-f| + l,|G|r2), r<R 


( 21 ) 


where cr = ±1 the upper sign corresponding to the A sublattice and the lower sign to the B one. 
We introduced the label <, > to separate between eignenergies of the two regions such that e< = FI 
and e> = E — Uq. Note that, the absolute value of G is introduced in order to consider the case of 
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positive and negative effective magnetic field. In addition, we note that the solution depend on 

the valley index r through (13). Such that </>+(r) and <^l(r) are envelope wave functions of A and B 
sublattices for the valley K, and 0^(r), 4>Z{r) for the valley K'. 

At this stage, it is natural to ask about the eigenvalues associated to the bound states of our 
system. To answer this inquiry we use the matching conditions of the spinors at the interface r = R 
to end up with the condition for the existence of bound states 

^ U{n>,\^-l\ + l,\G\R^) ^ ^ M{nU(-l\ + l,\G\R^) 

a- U{n>,\^ + l\ + l,\G\R^) f3- M {n<,\^ + l\ + 1,\G\R^) ^ ^ 

where the ratios a+ja- and l3+fl3- can be explicitly determined by making use of the coupled first- 
order differential equation (15). The sign of G depends on whether the oscillator frequency to is larger 
or smaller than the cyclotron frequency Uc- 

It is worthwhile investigating the basic features of of the above results and underlying their in¬ 
teresting properties. In this context, we would like to show how the valley degeneracy can be lifted 
by the presence of a constant magnetic field normal to the plane of the system together with a 2D 
Dirac-oscillator. This is of particular importance to form the valley filters, valves [30], or qubits [31] 
and spin qubits [32] in graphene. To this end, we consider three different cases corresponding to the 
relative strength of the magnetic field (cUc) with respect to the oscillator coupling (cj). 


3.2 Oscillating with same frequency 


We study the first case were the oscillator frequency is tuned to resonate with the cyclotron frequency, 
i.e. w ~ Wc, which it is equivalent to consider G = AA' with A' = ojc — oj. By requiring jA'j <C A we 
show that the hypergeometric functions reduce to bessel functions [34] 

U {n>,j,z) = (23) 

M {n<,j,z) = r(j) Jj-i (24) 

and the characteristic equation ( 22 ) takes the form according to ^ > 0 


VA2-(F;-t/)2^?-| 1 

( -{e- uyr) 



VE‘-X‘ 

(VA2 -{E- UfR) 




or ^ < 0 


x/x^-iE-uy^-i-l 1 

[xjx^ -{E- UYr) 

' E-U 

[yw^^Fdj 


yx^ -{E- UYr) 


[x/E'^ - A2 r) 


(26) 


We notice that there is a mapping between both of the two last equations, which is assured by the 
change ^ and r —)> —r. In general, these two equations can not be solved in the closed form. 

Results for the energy levels as a function of the dot radius are shown in Figure 2 for ^ | and 

U = 1.5 A. The black solid and blue dot dashed curves correspond to the K and K' valleys, respectively. 
Note that the case ^ for the same valley is also presented in both curves in Figure 2, which can 
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Figure 2: Energy levels of a circular quantum dot as a function of the dot radius i? for ^ ^ and 

U = 1.5 A, with A = hvp/X. The energy levels corresponding to the K and K' valleys are shown by the 
solid curve (r = 1) and the blue dot dashed curve (r = —1) , respectively. 

be discussed by taking into account of the symmetry r) = E{—^, —r) exhibited by our system. 
It may be noted that the solution for the two curves are different for the same valley. This can be 
explained by the fact that for the same valley we have E{^) ^ E{—^). In fact, we can easily see 
from (25) and (26) that the energy verify E{^,t) = E{—^,—t) and E{^,t) / E{—^,t), which helps 
us to study the valley degeneracy by a Dirac oscillator. However, the results show that the valley 
degeneracy is not broken, in this case, if we include both sign of Our results are in agreement with 
previous calculations of bound-state energies for zero magnetic field [35]. In summary, it is clearly 
seen that when the oscillator frequency is tuned to resonate with the cyclotron frequency, the Dirac 
oscillator annihilates the effect of the magnetic field on our system and we find ourselves in the case 
of zero effective magnetic field. At zero effective magnetic field, the degeneracy of the levels in the 
two valleys is clearly displayed. 

3.3 Strong magnetic field 

Now, we study the case when the cyclotron frequency is much larger than the oscillator frequency, i.e. 
Wc ^ w, then s = 1 and G ~ XtOc- It may be noted that the effect of the oscillator interaction is to 
produce an effective magnetic field along the positive 2 ;-direction that produces the effective quantized 
Landau levels. 

In this case we arrive, after applying the boundary conditions and finding both ratios, at the 
following characteristic equation for ^ > 0 

E-U+ tX U{n^+^ + l,( + l,GR‘^) _ ^ 

4G u{n> + GR^) ~ E - rX 

and for ^ < 0 

1 U{n>,-^ + lGR^) _ e-2 

E-U + tX U{n^ + !,-( +I,GR^) E + tX 


M (n^ + ^ ^+ ^, GR^) 

M{n< +^ +1,^ + 1,GR^) 


(27) 


M{n<,-^ + lGR^) 

M(n< + l,-e + |,G7?2) 


(28) 
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where the used quantum number is given by 


n 


<,> _ 

1 ~ 


- 4 .> 

4|G| • 


(29) 


Figure 3 shows the energies of the bound states in a graphene QD as function of the dot radius R 
for small effective magnetic field. In this Figure, we show the low-lying bound states for (^ = zb| 
with U = 105/17 = 205 (Figure 3 (a)/(b)) and A = 105. We can clearly show that for zero effective 
magnetic field the degeneracy of the levels in the two valleys is not broken. This is in agreement with 
the results found in the hrst case when the oscillator frequency is tuned to resonate with the cyclotron 
frequency. However, by increasing the effective magnetic field, the valleys degeneracy is broken. This 
allows us to conclude that the oscillator interaction, which allow to produce an effective magnetic 
held, can be used to control the orbital degeneracy. In addition, we show that the number of the 
bound states depend on the electrostatic potential. In fact, by increasing U the number of bound 
states decrease. 




Figure 3: Energy levels of a circular quantum dot as a function of the dot radius R with A = 10 5, 
[/ = 10 5/ [7 = 20 5 ((a)/(b)) and we have introduced the quantum dot level spacing 5 = hvp/R- The 
energy levels corresponding to the K and K' valleys are shown by the solid curve (r = 1) and the dot 
dashed curve (r = —1), respectively. We show the low-lying bound states for ^ ±|. 

In Figure 4 we show the energies of the bound states in a graphene QD as function of the dot 
radius R for large effective magnetic held for ^ ±|, ±|. The solid and dot dashed curve 

correspond respectively to the K and K' valleys. The results show that for large effective magnetic 
held the energy levels converge to the bulk Landau Levels. In addition, the states that are degenerate 
for zero effective magnetic held correspond to opposite values of the angular momentum ^ in different 
valleys. We note that our results are in agreement with previous calculations of bound-state energies 
in presence of magnetic held applied perpendicular to the graphene plane [35]. 





































RV^ 


Figure 4: Energy levels of a circular quantum dot as a function of the dot radius R with m = 10 <5, C/ = A 
and we have introduced the quantum dot level spacing 6 = hvp/R. The energy levels corresponding to the 
K and K' valleys are shown by the solid curve (r = 1) and the dot dashed curve (r = —1) respectively. 
We show the bound states for large effective magnetic field with ^ ±|, ±|. 


3.4 Weak magnetic field 


In the third case we suppose that the cyclotron frequency is much smaller than the oscillator frequency, 
namely tUc w, with G ~ —Xu, or s = — 1 . It may be noted that in this situation the effect of the 
oscillator interaction is to produce an effective magnetic field along the negative z-direction. In fact, 
this can be done by interchanging the confinement frequency u with the cyclotron one Uc- 

By using the above consideration, we obtain from (22) the characteristic equation for the allowed 
eigenenergies E of the QD. Indeed, for (^ > 0 we have 

1 U{n>,C + i-GR^) ^ M{n<,C+l-GR^) 

E-U-tX U{n>+ 1,^ + 1,-GR^) E - rX M{n< + 1,^ + ^,-GR"^) ^ ’ 


and for ^ < 0 reads as 

E-U-tX U{n^-C + l,-^ + l,-GR‘^) _ M(nf - ^ + i,-^ + 5 , 

4G u{n> -C + l,-^ + l,-GR^) ~ E + tX M(n< - ^ + i,-? + f, 

It is important to note that (30) is obtained from (28) by the following mappings 


-GR^) 

-GR^)' 


(31) 




-T, 


G 


-G 


(32) 


which can also be used to derive (31) from (27) in similar way as before. Moreover, the bound-state 
levels for negative effective magnetic field can be found by using the symmetry 


E{C,r,G) = Ei-^,-T,-G). (33) 

In both cases (positive and negative effective magnetic fields), the valley degeneracy is controllably 
broken by the presence of the effective magnetic field corresponding to our system. 
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4 Conclusion 


We considered the Dirac equation in (2+l)-dimensions in the presence of a constant magnetic field 
together with a two-dimensional Dirac-oscillator potential coupling. Indeed, using a similarity trans¬ 
formation, we formulated our problem in terms of the polar coordinate representation that allows us 
to handle easily the exact relationship between spinor components. Then, we got the solutions to the 
energy spectrum of the Dirac equation for a gate-tunable graphene QD. 

We discussed our results based on different physical parameters. In fact, we considered three cases, 
whether the cyclotron frequency is similar or larger/smaller compared to the oscillator frequency. This 
suggests defining an effective magnetic field that produces the effective quantized Landau levels. Our 
results were employed to discuss three important limiting cases those concern the weak, strong and 
fine tuned magnetic held situations. 

We studied the effective magnetic held dependence of energy levels. Indeed, we showed that the 
degeneracy of the valley is controllably broken by the effective magnetic held. In the hne tuned 
magnetic held cases (zero effective magnetic held), we found that the degeneracy of the levels in the 
two valleys is clearly displayed. Then, we analyzed in detail the case when the cyclotron frequency is 
much larger than the oscillator frequency (strong magnetic held). It has also been shown that at zero 
effective magnetic held, the valley degeneracy are not broken. However, by increasing the effective 
held we showed that the valley degeneracy are broken. 

In addition, when the magnetic held is turned off the problem becomes a pure Dirac oscillator, 
which leads to the creation of a negative effective magnetic held. These results showed that it is 
possible to control both spin and valley degeneracy by the effective magnetic held corresponding to 
our system. One can note that our results have an importance in forming valley hlters as well as spin 
qubits in graphene quantum dots. 
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